Quantum and classical criticalities in the frustrated two-leg 

Heisenberg ladder 



Mohamed Azzouqj and Brandon W. Ramakko 
Laurentian University, Department of Physics, 
Ramsey Lake Road, Sudbury, Ontario P3E 2C6, Canada. 
(Dated: August 21, 2007) 

This talk was about the Irustration-induced criticality in the antiferromagnetic Heisenberg model 
on the two-leg ladder with exchange interactions along the chains, rungs, and diagonals, and also 
about the effect of thermal fluctuations on this criticlity. The method used is the bond mean-field 
theory, which is based on the Jordan-Wigner transformation in dimensions higher than one. In this 
paper, we will summarize the main results presented in this talk, and report on new results about 
the couplings and temperature dependences of the spin susceptibility. 
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FIG. 1: The two- leg ladder showing the couplings along the chains, rungs, and diagonals is 
displayed. 

I. INTRODUCTION 

We use the bond- mean-field theory (BMFT), which is based on the Jordan- Wigner 
(JW) transformation, to study the quantum criticality phenomenon in the frustrated an- 
tiferromagnetic (AF) two-leg Heisenberg ladder, and the effect of temperature on this 
criticalityii^ii^iiii^. This method has been applied to the Heisenberg single chain, two-leg 
ladder, and three- leg ladder without frustration with excellent results^. When the diagonal 
interaction is varied the two-leg ladder system can undergo a quantum phase transition 
between two of three distinct non-magnetic quantum spin liquid states; the Neel-type 
(N-type) state, ferromagnetic-type rung (R-type) state, and ferromagnetic-type chain (F- 
type) state. These states are characterized by ferromagnetic spin arrangements along the 
diagonals, rungs, or chains respectively. The BMFT is a mean-field theory that is based 
on the spin bond parameters. The latter are related to the spin-spin correlation function 
{S~Sj'), with i and j labeling two adjacent sites in the direction where this correlation 
function is calculated. All quantities (5^ ), with a — x, y, z, are zero in BMFT, implying 
the absence of any sort of long-range magnetic order. 

The Hamiltonian for the spin-^ two- leg ladder with diagonal interactions is written as 

N 2 N N 

H = J • Si+ij + Jj^ • Si^2 + ' + • Si^2), (1) 

i j—l i i 

where J is the coupling along the chains, J± the coupling along the rungs, and Jx the 
coupling along the diagonals as seen in Fig. [TJ The index i labels the position of the spins 
along the two chains, each of which has TV sites, and j labels the chains. As usual, Sij is 
the spin operator. 

The frustrated two-leg ladder has been studied numerically using the Ising and dimer 
expansions^, the Lanczos diagonalization technique^'^'^'i^, and the density-matrix renor- 
nialization group (DMRG) ^°i^^'^^i^^'^'* . It has also been studied analytically using the 
bosonization2ii^, the valence-bond spin wave theory!^, a non-perturbative effective- field 
theoryii, the non-linear sigma model"., the reformulated weak-coupling field theoryi^, 
and the Lieb-Mattis theoreni2£. None of these analytical works addressed the issue of the 
phase diagram in all regimes including the weak, intermediate, and strong coupling lim- 
its. Within the BMFT, the phase diagram and its temperature dependence can be easily 
examined in all these regimes. 

This paper is organized as follows. In Sec.|TTl we explain how the BMFT, which is based 
on the JW transformation in dimensions higher than one, is applied to our Hamiltonian. 
The Quantum and classical critical behaviours, and the spin susceptibility are discussed 
in Sec. IIIII In Sec. IIVI conclusions are reported. 



3 



II. METHOD 

The JW transformation for the two-leg Heisenberg ladder is defined as^ 

^^J^ '^ ' ^iJ^^J^ 
i-1 2 

= 7''[^^ndj] for chain 1, 

d=Qf=l 
i-l 2 

0i,2 = 7r[^ ^ rirfj' + rii^i] for chain 2. (2) 

The cjj- operator creates a spinless fermion at site while annihilates one, and 

Uij is the occupation number operator at that site. The phases (f>ij are chosen so that at 
the the spin operators commutation relations are preserved. 

After applying the JW transformation ([2]) to the Hamiltonian ([T]) we get 

N N 

^ =2 E(4ie^""-'c,+i,i + cl ^e— '+^.^c.+i,2 + H.c.) + --f + H.c.) 

i i 
N N 2 

+ ^E(4ie"'""^'"""+^'^^c,+i,2 + 4i,ic,2 + H.C.) + jJ2Y.(''u - 2^("^+^'^' ~ 2^ 

i i j=l 

N ^ , N 



+ ("»+!,!- ^)(^i,2 - ^)]. (3) 

In BMFT, the interacting terms of the JW fermions are decoupled using the spin bond 
parameters. This approximation neglects fluctuations around the mean field points; (O — 
{0)){0' — {O')) « 0, where O and O' are any operators which are quadratic in and 
This yields 

00' ^ {0)0' + 0{0') - {0){0'). (4) 

To apply BMFT we introduce three mean-field bond parameters; Q in the longitudinal 
direction, P in the transverse direction, and P' along the diagonal. These can be inter- 
preted as effective hopping energies for the JW fermions^ in the longitudinal, transverse 
and diagonal directions, respectively: 

Q ^ (cjjcj+i^j), P ^ P' = (c,+ijc^^i). (5) 

Keeping in mind that there is no long-range order-1 so that (S^j) = (c|-^Ci.i) ^ 1/2 = 0, 
the Ising quartic terms in equation ([3]) can be decoupled using the Hartree-Fock approxi- 
mation ([5]), and the bond parameters ([5]) as follows: 



4,1(^^,1- 2) 2) "S^iC^+i,! + "5*4+1, i*^'.! + l'?!^ (6) 

for the Ising interaction along the chain 1. Similar equations can be obtained for the 
interactions along the other directions. Next, we write the Hamiltonian ([3]) using the 
three different spin configurations in the right panel of Fig. [21 These configurations are 
instantaneous (not static) configurations in which adjacent spins in any direction keep on 
average the same relative orientations with respect to each other, but fluctuate globally 
on a time scale determined by the strongest coupling constant so that any kind of long- 
range magnetic order is absent. These fluctuations are a consequence of the spin quantum 
fluctuations. The three competing configurations of Fig. [2]lead to three different quantum 
gapped spin liquid states, each characterized by its own short-range spin correlations and 
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FIG. 2: In the left panel, the three possible ground states of the system in the Ising limit, namely 
the Neel state, the ferromagnetic chain state, and the ferromagnetic rung state are drawn. In the 
right panel, the labeling of sublattices corresponding to the short-range spin orders that replace 
the long-range ones are shown for the Heisenberg limit. 



symmetry. We also choose to place an alternating phase of tt along the chains so that 
the phase per plaquette is tt—. This phase configuration is used to get rid of the phase 
terms in the Hamiltonian. We also set Qij = Qe'*'-^ where Q is site independent^. Here 
is the phase of the bond along the chain; i.e., if = on a given bond, then it is 
equal to tt on the adjacent ones. This is necessary in order to recover the proper result in 
the limit Jx and J± becoming zero, in which we get an energy spectrum comparable to 
that of des Cloiseaux and Pearson^^ for the spin excitation in a single Heisenberg chain, 
E{k) = f J|sinfc|. 

For each state the Hamiltonian is written in the Nambu formalism and the matrix is 
diagonalized to obtain four eigenenergies (for each of the states). The details can be found 
in Refi^. The eigenenergies for the N, F, and R-type states are given respectively by 



EN{k) = ±Jx 1 cos fc ± \ Jl sixC k 
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J2 

Epik) = ±Ji cosfc ± y sin^ k + 



ER{k) ^ ±^ ± y^j2sin2/c + j2^cos2fc, (7) 

where Ji = J(l + 2Q), J_ii = J_l(1 + 2P), and Jxi = Jx (1 + 2P')- The free energy per 
site is 

= + ^ + ./xP'^ - ^ E E + ^""''''^'^]' (8) 

k p=l 

where Ep{k) is one of the four eigenenergies in any state. The minimization of the free 
energy with respect to the bond parameters leads to the following set of self-consistent 
equations: 

k p=l 



4:NJ± dP 

k p—1 



P' = ^—yy^^Mnp[EJk)], (9) 

fc p=i 



which are solved numerically, except in the high-temperature regime, where analytical 
results are obtained. 
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FIG. 3; The zero-T (a) and finite-T (b) (ai, a2)-phase diagrams are shown. The zero-T one is 
compared with the Lanczos- method'^ and the Ising hmit. The boundaries are between the N-type 
state (N), R-type state (R), and the F-type state (F). 




FIG. 4: The uniform spin susceptibility xiT) is plotted as a function of temperature for h = 
for several coupling constants. 



III. RESULTS 



The free (ground-state) energies of all three states are calculated as functions of the 
coupling constants and compared. From thermodynamic considerations the state with the 
lowest free energy is the stable one, and whenever free energies cross a phase transition 
takes place. Since only the ratios of the couplings are important we define ai = J±/ J and 
ct2 = Jx/ J- In this way we have obtained the zero and finite temperature phase diagrams 
which can be seen in Fig. [3l The agreement between the Lanczos method data and our 
results is very good, a fact that indicates that the present mean-field treatment is accept- 
able. The line at q;2 = 1 is exact and its placement is a consequence of the Hamiltonian 
symmetry with respect to exchanging J and Jx . The quantum phase transitions (zero 
temperature) found here using BMFT are first-order ones for all values of a2- As tem- 
perature increases the R-type state decreases in size. The sizes of the N-type and F-type 
phases increase with temperature. For any set of coupling values in the shaded region of 
Fig. a classical (thermally induced) first order transition from the R-type state to one 
of the other states occurs-. 

The uniform magnetic susceptibility x(T) is calculated following the same method as 
Refi^. It is plotted as a function of temperature in Fig. |4] . The Heisenberg model on a 
chain (ai = 0, Q!2 — 0) is gapless whereas when ai and/or q;2 7^ an energy gap 
opens up in the low-energy excitation. This is why at zero temperature x 7^ only for the 
chain, and an exponentially decreasing susceptibility is obtained for nonzero ai or a2- For 
ai = 0.6 and a2 — 0.5 there is a sudden transition from the R-type state to the N-type 
state at kgT/ J = 0.39. For example, for ai = 1.25 and 02 = 2 there is a sudden transition 
from the R-type state to the N-type state at fc^T/J — 0.64. This sudden transition results 
in a discontinuity in x- These findings remain to be confirmed by others means. 
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IV. CONCLUSION 



In this talk, quantum and classical critical behaviours in the frustrated antifcrromag- 
netic two-leg ladder were presented. The method of calculation, which is based on the 
Jordan- Wigner transformation and the bond-mean-field theory was explained. The zero- 
temperature phase diagram of this system was explained. It exhibits three quantum 
phases, characterized all by an energy gap and absence of magnetic order. These states 
are labeled Neel-type, Rung-type and Ferromagnetic-type chain states. Our zero-T results 
agree well with existing numerical data. When temperature increases for some sets of cou- 
pling values, the system undergoes a phase transition from the R-type state to the N or 
F-type state at a finite temperature. The finite temperature phase diagram was explained 
as well. In it, the size of the R-type state becomes smaller while the F-type state and the 
N-type state increase in size with increasing temperature. Our theory predicts a disconti- 
nuity in the spin susceptibility at a finite temperature for the sets of couplings where the 
finite-T transition occurs. 
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